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Much of our knowledge of galaxies comes from analysing the radi- 
ation emitted by their stars. It depends on the stellar initial mass 
function (IMF) describing the distribution of stellar masses when 
the population formed. Consequently knowledge of the IMF is 
critical to virtually every aspect of galaxy evolution. More than 
half a century after the first IMF determination^, no consensus has 
emerged on whether it is universal in different galaxies®. Previ- 
ous studies indicated that the IMF and the dark matter fraction 
in galaxy centres cannot be both universal^], but they could not 
break the degeneracy between the two effects. Only recently in- 
dications were found that massive elliptical galaxies may not have 
the same IMF as our Milky Way 8 . Here we report unambiguous 
evidence for a strong systematic variation of the IMF in early-type 
galaxies as a function of their stellar mass-to-light ratio, producing 
differences up to a factor of three in mass. This was inferred from 
detailed dynamical models of the two-dimensional stellar kine- 
matics for the large ATLAS 3D representative sample^ of nearby 
early-type galaxies spanning two orders of magnitude in stellar 
mass. Our finding indicates that the IMF depends intimately on 
a galaxy's formation history. 

As part of the ATLAS 3D projecP, we obtained integral-field maps 
of stellar kinematics for a volume-limited sample of 260 early-type (el- 
liptical and lenticular) galaxies. They were selected to be closer than 
42 Mpc and to have a its -band total magnitude Mr < —21.5 mag 
(M > 6 x 10 9 M ), as determined from the Two Micron All Sky Sur- 
vey (2MASS) at our adopted distances. Homogeneous imaging for all 
the galaxies in the r-band was obtained in major part from the Sloan 
Digital Sky Survey (SDSS) DR8 and completed with our own photom- 
etry. D 

For all galaxies, we constructed six sets of dynamical modelP, 
which include an axisymmetric stellar component and a spherical dark 



1 Sub-department of Astrophysics, Department of Physics, University of Oxford, Denys 
Wilkinson Building, Keble Road, Oxford 0X1 3RH 2 Gemini Observatory, Northern Operations Cen- 
tre, 670 N. A'ohoku Place, Hilo, HI 96720, USA 3 Department of Astronomy, Campbell Hall, Uni- 
versity of California, Berkeley, CA 94720, USA J Observatoire de Paris, LERMA and CNRS, 61 Av. 
de I'Observatoire, F-75014 Paris, France 5 Laboratoire AIM Paris-Saclay, CEA/IRFU/SAp CNRS 
Universite Paris Diderot, 91191 Gif-sur-Yvette Cedex, France 6 Department of Astrophysics, Univer- 
sity of Massachusetts, 710 North Pleasant Street, Amherst, MA 01003, USA 7 European Southern 
Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany 8 Sterrewacht Leiden, Leiden 
University, Postbus 9513, 2300 RA Leiden, the Netherlands 9 Universite Lyon 1, Observatoire de 
Lyon, Centre de Recherche Astrophysique de Lyon and Ecole Normale Superieure de Lyon, 9 av- 
enue Charles Andre, F-69230 Saint-Genis Laval, France ln Max-Planck Institut fur extraterrestrische 
Physik, PO Box 1312, D-85478 Garching, Germany ^Netherlands Institute for Radio Astronomy 
(ASTRON), Postbus 2, 7990 AA Dwingeloo, The Netherlands 12 Kapteyn Astronomical Institute, 
University of Groningen, Postbus 800, 9700 AV Groningen, The Netherlands 13 Max-Planck Institut 
fur Astrophysik, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany 14 Centre for Astrophysics 
Research, University of Hertfordshire, Hatfield, Herts AL1 9AB, UK 15 Centre for Astrophysics & 
Supercomputing, Swinburne University of Technology, PO Box 218, Hawthorn, VIC 3122, Aus- 
tralia 16 Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St. George Street, 
Toronto, ON M5S 3H4, Canada 17 Physics Department, New Mexico Institute of Mining and Tech- 
nology, Socorro, NM 87801, USA 



halo, and fit the details of both the projected stellar distributioiP and 
the two-dimensional stellar kinematics-^ (Fig. 1). While the shape of 
the stellar component can be inferred directly from the galaxy im- 
ages, the dark halo shape must be a free parameter of the models. 
We explored with the models a variety of plausible assumptions for 
the halo to test how these can affect our result. Our halo models in- 
clude as limiting cases a maximum-ignorance one, where the halo pa- 
rameters are directly fitted to the stellar kinematics, and some com- 
pletely fixe d ones, where the halo follows the predictions of numerical 
simulation J H l 13 l l4 [ A detailed description of the model parameters is 
provided in Table 1. The key parameter we extract from all the models 
is the ratio (M / X) sta rs between the luminosity (in the r-band) and the 
mass of the stellar component. As illustrated in Fig. 1, the availability 
of integral-field data is the key to accurately separate the stellar mass 
from the possible dark matter using dynamical models and determine 
(M/L) s t ars . In fact changes in (M/L) star s at the level expected for 
IMF variations cause dramatic changes to the quality of the model fits. 

We also measured the (M/L) pop of the stellar population by 
fitting^ the observed spectra using a linear combination of single stel- 
lar population synthetic spectrPof different ages (t) and metallicities 
([M/H]), adopting for reference a Salpeter-Q IMF (£(m) oc m x = 
m~ 2,3 ). The models adopt standard lower and upper mass cut-offs for 
the IMFofO.l and 100 Mq, respectively. We used linear regularization 
to reduce noise and produce smooth M(t, [M/H]) solutions consistent 
with the observations. The resulting (M/L) pop is that of the compos- 
ite stellar population, and excludes the gas lost during stellar evolution. 
If all this gas was retained in the galaxies in gaseous form, it would 
systematically increase (M/L) pop by about 30"%P^. However most of 
it is likely recycled into stars or expelled to larger radii. Although the 
results are cleaner using our full spectrum fitting approach^! similar 
conclusions are reached when the galaxies are approximated as one 
single stellar popula tion, o r when (M/L) pop is computed using differ- 
ent population codes^ 7 ^ and with a more traditional approach which 
only focuses on the strength of a few stellar absorption spectral lines. 
Systematic offsets of about 10% in (M/L) pop exist between the pre- 
dictions of different population models, for an identical set of assumed 
population parameters, with the adopted one being in the middle of the 
others. This sets the uncertainty in the absolute normalization of our 
plots. The random errors for our adopted population codJ 7 ^ were es- 
timated by applying the same spectral fitting approach to our integral- 
field spectroscopy data and to independent spectra obtained for a subset 
of 57 of our galaxies by the SDSS survey. We inferred an rms scatter 
of 12% in each individual (M/L) pop determination. 

The ratio between the dynamically-derived (M/X) s tars and the 
population-derived (M/L)salp, using a fixed Salpeter IMF, is shown in 
Fig. 2 as a function of (Mj L) 8 t ar8 . We compare the observed ratio with 
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Figure 1 

V rIns = \/V 2 + ct 2 for five galaxies representing a variety of shapes of the kinematics fields, and spanning a range in (MIL) 



Disentangling the stellar and dark mass with integral-field stellar kinematics. The top panels show the symmetrized SAURON stellar kinematics 

stars values. Here V is the mean 



stellar velocity and a is the stellar velocity dispersion. The middle panel is the best-fitting dynamical modePSwith a standarcP-Q dark halo (model b in Table 1). The 
bottom panel is a dynamical model where the (M / L) ata rs was fixed to be 0.65 times the best-fitting one. Where this decrease in (M / L) ata rs represents the change 
in mass between a Salpeter and Kroupa IMF. The other three model parameters, the galaxy inclination i, the orbital anisotropy /3 Z and the halo total mass M200. were 
optimized to fit the data, but cannot provide an acceptable description of the observations. This plots shows that, for a standard halo profile, the data tightly constrain 
both the dark matter fraction and (M / -L) s tars- The effect would be even more dramatic if we had assumed a more shallow inner halo profile. The contours show the 
observed and modelled surface brightness respectively. The values of (M/ L) Bta ,rs and the fraction of dark matter within a sphere with radius equal to the projected 
half-light radius R e are printed next to each panel. The galaxy names are given at the top. 
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Figure 2 | The systematic variation of the IMF in early-type galaxies. The six panels show the ratio between the (M/L) ata rs of the stellar component, determined 
using dynamical models, and the (M/L)s a i p of the stellar population, measured via stellar population models with a Salpeter IMF, as a function of (M/L) B tars- 
The black solid line is a loess smoothed version of the data. Colours indicate the galaxies' stellar velocity dispersion <j e , which is related to the galaxy mass. The 
horizontal lines indicate the expected values for the ratio if the galaxy had (i) a Chabrier IMF (red dash-dotted line); (ii) a Kroupa IMF (green dashed line); (iii) a 
Salpeter IMF (x = —2.3, solid magenta line) and two additional power-law IMFs with (iv) x = — 2.8and(v)n: = — 1.5 respectively (blue dotted line). Different 
panels correspond to different assumptions for the dark matter halos employed in the dynamical models as written in the black titles. Details about the six sets of models 
are given in Table 1. A clear curved relation is visible in all panels. Panels a, b and e look quite similar, as for all of them the dark matter contributes only a small 
fraction (zero in a and a median of 12% in b and e) of the total mass inside a sphere with the projected size of the region where we have kinematics (about one projected 
half-light radius R e ). Panel f with a fixed contracted halo, still shows the same IMF variation, but it is almost systematically lower by 35% in (M / L) s tars reflecting 
the increase in dark matter fraction. The two black thick ellipses plotted on top of the smooth relation in panel d show the representative lcr error for one measurement 
at the given locations. We excluded from the plot the galaxies with very young stellar population (selected as having H/3 > 2.3 A absorption). These galaxies have 
strong radial gradients in their population, which break our assumption of spatially constant M j L and makes both (M/ L)s a i P and (M/L) B tars inaccurate. 
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Tabic 1 | The axisymmetric dynamical models 



Panel in Fig. 2 Description of the model Fitted model parameters 

a Galaxy model in which the total mass traces the observed galaxy i, j3 z , (M / L) tota \ 

light distribution. Any dark matter, if present, follows the stellar 
distribution. 

b Galaxy stellar component embedded in a spherical standard dark i,P z , {M / L) Bta , IB , M200 

matter with inner density p(r) oc r 1 for radii r <C r s and 

outer density p{r) oc r~ 3 for r>r s . The halo total mass M200 is 

fitted, while r a is uniquely specified^ by A/200, 
c Model with a standard 1 1 halo contracted^ according to the observed i, f3 z , (Mj L) Bta , TB , M200 

galaxy stellar density. The halo mass is fitted, while r s is specifie (P 

by A/ 2 oo- 

d Model with a general halo inner density p(r) oc r 1 with fitted slope i, /3 Z , (M/L) s tars , 7, AZ200 

(—1.6 < 7 < 0) and total mass. The outer density becomes p(r) oc 

r~ 3 as in the standard hakP^at radii r 3> r s — 20 kpc. 
e Model with a fixed standard hakP with M 2 oo specified^! by the i,/3 z , (AZ/L) stars 

measured galaxy stellar mass and r s specified*^ by Af2oo- 
f Model with a fixed standard halrP contracted^ according to the i,P z , (M / L) Bt ^ rB 

observed galaxy stellar density. A/200 is specifie cPby the measured 

galaxy stellar mass and r s is specified^ by A/200 ■ 



the expected one if the galaxy had the 'light' KroupiP or ChabrieP 
IMFs, which are similarly deficient in low mass stars; the 'standard' 
Salpeter IMF, which is described by a simple power-law in stellar mass 
with exponent x — —2.3; and two additional 'heavy' power-law IMFs 
with x — —2.8 and x — —1.5 respectively. The last two IMFs predict 
the same (M / L) pop . But while in the first case the stellar population 
is dominated by dwarf stars, in the second case the large (AI/L) pop is 
due to stellar remnants: black holes, neutron stars and white dwarfs. 
The dynamical mass measurements do not constrain the shape of the 
IMF directly, but only the overall mass normalization, and for this rea- 
son cannot distinguish between the two cases. 

The results from all sets of dynamical models are consistent with a 
similar systematic variation of the IMF normalization, by up to a fac- 
tor of three in mass. A clear trend is visible in particular for the most 
general of our set of models (panel d of Fig. 2), which makes virtu- 
ally no assumptions on the halo shape but fits it directly to the data. 
However similar trends are visible for all our plausible assumptions for 
the dark halo mass and profile as predicted by numerical simulations. 
This shows that, although our result does not depend on the correct- 
ness of the assumed halo model, it is entirely consistent with standard 
model predictions for the halo. For increasing (A//L) sta rs the normal- 
ization of the inferred IMF varies from the one of Kroupa/Chabrier up 
to an IMF more massive than Salpeter. The trend in IMF is still clearly 
visible when selecting a subset of 60 galaxies lying outside the Virgo 
galaxy cluster, with the most accurate distances and the best models 
fits. This shows that it cannot be due to biases in the models, distances, 
or to effects related to the cluster environment. The knee in the rela- 
tion at (M/L) B tars ~ 6 (r-band) shows that the lowest (M / L) s t a i a 
values mainly reflect the age and metallicity of the population (with 
younger ages or lower [M/H] decreasing (A//L) pop ), while the largest 
(Af/L) at ars values mainly reflect their dwarf or remnants dominated 
IMF. The models with contracted halos show the same IMF trend as the 
other models. However contracted halos predict too little stellar mass 
for many of the galaxies with low (A//Z/) st ars, even for the 'lightest' 
Kroupa/Chabrier IMF. This suggests that contraction may not happen 
in most real galaxies, in agreement with recent numerical simulations^. 

Our result reconciles a number of apparently contradictory results 
on the normalization of the IMF accumulated in the past decade. The 
Rroupa/Chabrier-like normalization at low (Af/L) s tars agrees with the 
one inferred for spiral galaxies^. The ATLAS 3D project discovered 



that early-type galaxies with the lowest (M/L) pop resemble spiral 
galaxies with their gas and dust removed^ and thus a similarity of 
IMF should be expected. The Kroupa/Chabrier normalization is also 
consistent with previous findings that this normalization i s requi red 
not to over-predict the mass of early-type galaxies as a clasJ™ A 
Salpeter normalization at larger (M/Z/) st ars is consistent on average 
with results from strong gravitational lensing^, which are restricted to 
the galaxies with the largest velocity dispersions (<r > 200 km s _1 ). 
Finally, the normalization more massive than Salpeter for some of the 
galaxies with large (A//L) s tars is broadly consistent with the finding 
from the depth of spectral features of eight massive galaxiePthat indi- 
cate they must be dominated by a population of dwarf stars. 

If instead the largest (M/L) pop were due to stellar remnants, our 
results would be consistent with indirect arguments based on the rela- 
tion between the colour of a stellar population and its fraction of ioniz- 
ing photons, suggesting an I MF slo pe becoming flatter for more mas- 
sive and star forming galaxieP^H. However our result is difficult to 
compare directly, due to the large difference in the sample selections. 
Moreover these studies— 1 — 1 measure the instantaneous IMF, when the 
stars are forming, while all previous studies we mentioned, including 
the one in this Letter, measure the 'integrated' galaxy IMF (IGIMF) re- 
sulting from the cumulative history of star formation^ and evolutionary 
mechanisms that the galaxy has experienced. 

The discovered trend in IMF is also consistent with previous find- 
ings that the total M/L in the centre of galaxies varies by at least a 
factor of two more than one would expect for a stellar population with 
constant dark matter fraction and a universal iMF=4 But various previ- 
ous attempts could not distinguish whether t he mas s discrepancy was 
due to non-universality of dark matter or iM I* 4 * 7 ! 30 ! The studies were 
limited either by small samples or non optimal data^ or employed 
simplified g alaxy m odels that could bias the quantitative interpretation 
of the resultP^M! \y e finally resolve both of these limitations. 

Our study convincingly demonstrates that the assumption of a 
universal IMF, which is still adopted in nearly every aspect of galactic 
astrophysics, stellar populations and cosmology, is inconsistent with 
real galaxies. The results presented pose an interesting challenge to 
galaxy formation models, which will have to explain how stars 'know' 
what kind of galaxy they will end up inside. A possible explanation 
would be for the IMF to depends on the prevailing physical conditions 
when the galaxy formed the bulk of their stars. Although galaxies 
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merge hierarchically, there is growing evidence that present-day 
massive early-type galaxies formed most of their stars in more intense 
starbursts and at higher redshift than spiral galaxies. This could lead to 
the observed difference in IMF. Unfortunately there is no consensus 
in the theoretical models for how the IMF should vary with physical 
conditions. A new generation of theoretical and observational studies 
will have to provide insight into which physical mechanisms are 
responsible for the systematic IMF variation we find. 
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